Event shapes for hadronic final state: experimental review 
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Abstract. Analysis is presented for first moments of collective observable distributions in two-jet events 
for various interaction types and for wide initial energy range. These observables include sphericity, thrust, 
components of transverse particle momentum, alignment and planarity. Database of experimental results 
created in the framework of the paper includes data for all interactions. Energy dependencies of average 
values for collective observables excepting of transverse momentum components show universal behavior for 
various interactions. Thus energy dependence for these parameters supposes the descriptions by functions 
which are universal for different interaction types. Particle transverse momentum as well as it's components 
increase faster for e + e~ interaction with growth of y/s, than that for other interactions. Empirical analytical 
functions are suggested for description of energy dependence for all collective observables under study with 
exception of infrared-stable thrust variable. Energy dependence for average thrust is compared with QCD 
predictions included perturbative part and analytical phenomenological corrections which account for non- 
perturbaive effects. Dispersive model and single dressed gluon approximation are considered for description 
of energy dependence of first moment of thrust distribution and estimation of strong coupling constant for 
various interactions as well as for joint sample. The dispersive model allow to describe average thrust vs 
initial energy in wide range of y/s down to strongly non-perturbative domain y/s ~ 2 — 3 GeV at qualitative 
level at least. Study of event shape observables allows to obtain estimations of as{Mz) which are in 
reasonable agreement both with world average value and with results extracted in the framework of other 
methods. Using suggested analytic approximation functions estimations of values of collective parameters 
under study have been obtained for present and future facilities. In TeV energy domain average values of 
collective observables either depend on y/s weakly or do not depend on initial energy at all within errors. 
Thus, the TeV scale can be considered as an estimation of the low bound of asymptotic region for traditional 
collective parameters. Usually, multiplicity dependence of collective observables under considered agree 
with power function in energy domain y/s < 12 GeV at qualitative level at least. Behavior of sphericity 
vs multiplicity and comparison of experimental results with model calculations allow to suggest that the 
universal estimation of the low energy boundary for experimental appearance of event jet structure in 
multiparticle production processes is y/s~2 ~ 3 GeV. 

PACS. 13.87.-a Jets in large-Q 2 scattering - 12.38.-t Quantum chromodynamics 

1 Introduction into the fundamental workings of QCD, since they con- 

tain within themselves the signatures of QCD at both 

T . . , , . , weak and strong coupling. Jets therefore test our under- 

m various interactions and m wide initial energy do- ,. r ^^5. . , r . f1l T , 
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.,, . . . . i m? i i • cays ot massive particles like top quarks. In order to map 
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titative definition ot hadronic jet tor various experiments 

Since we are not ab e to directly measure colored ob- ig & difficult problem To a first approximation at en ough 

jects (quarks or gluons), but only hadrons and their de- high energieSj a jet can be thought of as a hard parton 

cay products, a central issue for every experimental test of th&t hag undergone soft and co llinear showering and then 

quantum chromodynamics (QCD) is establishing a corre- hadronization . Good jet definitions have to be infrared 

spondence between observables obtained at the partomc and coUinear gafej simple tQ uge in theoretical and experi _ 

and the hadronic level. Hadronic jets provide a window mental contextS; applicable to any type of inputs (parton 
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ergy deposits in the detectors) and lead to jets that are 
not too sensitive to non-perturbative effects [2]. The re- 
view of investigations of jet definitions is given in [3] (for 
e+e~ annihilation) and in [4, 5] (for hadronic - pp / pp - 
collisions). 

On the one hand, one can study observables that de- 
pend on specifying the actual number of jets in the final 
state by defining a jet algorithm, which sets criteria for 
what constitutes a jet. There are two main classes of such 
algoritms: cone algorithms, extensively used at hadron 
colliders, and sequential recombination algorithms, more 
widespread at e + e~ and ep colliders. On the other hand, 
one can extract much useful information about the struc- 
ture of the final state from simpler observables called event 
shapes or collective variables, which do not depend on a 
jet algorithm, but are simple functions of the momenta of 
all final state particles included in jet analysis. 

Event shapes are used for many purposes. Some col- 
lective characteristics can be calculated perturbatively in 
terms of QCD colored objects and compared with hadronic 
final state measurements at high energies [6]. The LEP 
experiments and their high quality data have prompted 
much theoretical progress in the understanding of QCD 
radiation and in the development of appropriate tools. 
This includes Sudakov resummation, renormalon resum- 
mation and parametrization of power corrections [7]. At 
the same time, experimental uncertainties and non-per- 
turbative effects are small in the regime of high trans- 
verse momenta, which makes possible precision tests of 
perturbative QCD and the extraction of the strong cou- 
pling constant. On the other hand, collective observables 
are directly sensitive to the hadronization process, provid- 
ing an opportunity to learn more about the confinement 
mechanism from experiment [7]. Understanding the na- 
ture of hadronisation is a challenge, which is crucial im- 
portant for all high energy experiments involving hadrons 
[8]. Thus the study of collective observables in hadronic 
collisions and for wide energy domain seems the impor- 
tant task and promising tool for investigation of problem 
of confinement, hadronization mechanism and correspond- 
ing analytical models. At the same time that hadronic jets 
provide insights into QCD, they are key elements in signa- 
tures for new physics beyond the Standard Model. Some 
traditional event shape observables can be used for search 
of such signatures. The detail investigations of the spheric- 
ity and thrust distributions are promising tools for sharp 
and robust discrimination between the SUSY and some 
models with large extra dimensions which should be util- 
isable either at ILC or at CLIC wherever these scenarios 
become accessible [9] . At a hadronic colliders such as the 
LHC, most interesting processes involve final states with 
jets of hadrons. The data available to define and recon- 
struct the jets are energy depositions in the calorimeter 
and charged particle tracks. Ideally one can use this ex- 
perimental data to identify an underlying hard scatter- 
ing event that involved QCD partons in the weakly cou- 
pled high energy regimes or, alternatively, heavy standard 
model particles such as top quarks or newly produced non- 
standard model particles. Distinguishing the former from 



the latter is of course crucial to our ability to find and 
analyze new physics at the LHC [10]. Thus traditional 
collective variables can be useful for search of some ex- 
otic events in hadronic collisions at (ultra)high energies. 
As expected ordinary QCD events can be distinguished 
from events involved decays of new particles because val- 
ues of some event shapes will be higher (or closer to the 
spherical limit) for latter case than that in the former one. 
In particular, higher sphericity of an event in relation to 
any other event of interactions described by the Standard 
Model is of the general regularities for multidimensional 
black hole production at accelerators [11]. But one need 
to emphasize that the challenge of calculating QCD back- 
ground events in the hadron collider environment of the 
RHIC, Tevatron or the LHC remains formidable so far. 

This paper focuses at experimental study of event sha- 
pes namely. At present the pp colliders RHIC and LHC 
as well as Tevatron allows to study the geometry of final 
state in new energy ranges at high precisions. Taking into 
account importance of study of hadronic interactions but 
not only collisions with lepton beams this paper pais at- 
tention, especially, to event shapes in hadron-hadron and 
hadron-nuclear reactions. 

The paper is organized as follows. In Sec. 2, definitions 
of collective variables under study are described. Database 
created in the framework of the paper is shown and dis- 
cussed in Sec. 3. The Sec. 4 devotes to discussion of energy 
dependence of collective variables, comparison it with var- 
ious approximations and extracted estimations of strong 
coupling constant. Section 5 demonstrates a multiplicity 
dependence of jet characteristics and in Sec. 6 some final 
remarks are presented. 

2 Observables in analysis of hadronic final 
state properties 

Event shape variables are functions of the four mo- 
menta in the hadronic final state that characterize the 
topology of an event's particle / energy flow. They are 
sensitive to QCD radiation (and correspondingly to the 
strong coupling) insofar as gluon emission changes the 
shape of the particle / energy flow [2]. The relationship 
of pattern of collective particle / energy flow with mea- 
sured parameters (momenta) of final state particles can 
be described by momentum tensor A a @ which is defined 
generally in the following way: 

N 

A afi = Y^WiPiPf), a,/3 = x,y,z. (1) 

Here N is the number of particles separated for collective 
analysis in event, Wi - some weight of i^ 1 particle. Both 
the collective flow of particles and energy can be investi- 
gated at choosing the suitable Wi. One need to emphasize 
that the form of the momentum tensor varies for different 
papers (see, for example, [12 16]). 

One of the most important and global properties of 
hadrons produced in various collisions in wide energy range 
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is their jet nature which can be characterized by the some 
collective observable. The sphericity (S) and thrust (T) 
are the most popular and widely used variables for study 
of jet behavior of the secondary particles. The sphericity 
is defined as the following [12, 13]: 



N 

3 

S = min S (n) = - min ^ 

2 " ^ 2 



(2) 



Here p is the transverse momentum of the i ^ particle 
with respect to the sphericity axis, which is chosen in such 

N 

a way that the sum ^2 Pj_i w ith respect to it is minimal; 

i=i 

Pi is the momentum of the particles in the center-of-mass 
system, and the sum extends over all particles separated 
for jet analysis. Since sphericity is quadratic in the mo- 
menta it cannot be calculated in perturbation part of the 
QCD, however it has remained a popular jet measure be- 
cause it can be easily determined [17]. A measure of the 
event structure that uses the linear momenta is the thrust 
[18 - 20] 

N 

E Ipiiil 



T = maxT(n) = max ■ 



JV 

E\Pi 



(3) 



N 



The thrust axis being chosen to maximize ^2 \p\u\, the 

i=i 

sum of the components of the momenta parallel to the 
thrust axis of all particles in an event selected for jet anal- 
ysis. 

Sphericity tensor S a ^ can be derived as special case of 
general form ([T]) at Vi : Wi = const. Usually constant is 

N 

chosen equal to unit or as follow const -1 = ^2 pf. The 

i=l 

present work is used the first case of constant and the 
matrix of sphericity tensor is the following: 



N 



x, y, z. 



(4) 



Because the sphericity tensor uses the momenta of par- 
ticles quadratically, the high-momentum particles in an 
event will contribute more strongly to observables derived 
from this tensor than to those which use the momenta lin- 
early [21]. To eliminate false effects which may be caused 
by leading particles, the reduced-momentum matrix [22] 
can be also used to study the configurations of inelastic 
collisions. 

Study of shape of hadronic final state in various in- 
teractions is carried out in the system of coordinates of 
the "principle axes" of the event, which are the symmetry 
axes for a given configuration of secondary-particle mo- 
mentum vectors (PigH]). This coordinate system has some 
special advantages for a more detailed study of the con- 
figuration of events in 3-momentum space. Conversation 



of this coordinate system was accomplished by standard 
procedure of diagonalization of a matrix (j4]) . Let n% , ri2 , 
and «3 be the unit eigenvectors of the tensor @ associ- 
ated with the eigenvalues A,, , which are ordered such that 

3 

Ai > A2 > A3 and A^ = J2 ^i- F° r this choice the vector 

i=l 

ni determines the direction of the greatest extent of the 
event in momentum space; 712 determines the direction 
of greatest extent in the plane perpendicular to ni, and 
ns determines the direction of the greatest compression in 
this plane. The axis coincided with direction of eigenvector 
ni is jet axis in the assumption that the certain event show 
two-jet structure. The diagonal elements of the matrix Q 
represent the sum of the squares of the components of the 
secondary-particle momenta in the coordinate system of 
the principal axes. If the diagonal elements are related as 
Ai ~ A2 ~ A3, then the event is spherically symmetric; if 
Ai ^> A2, Ai S> A3, and A2 ~ A3, the event has cyllindrical 
symmetry; if Ai ~ A2 and Ai, A2 3> A3, then the event 
is planar, recalling a disk in shape [23]. Forms of events 
at various values of traditional collective observables are 
presented, for example, in [20]. 




Fig. 1. System of coordinates of the "principle axes" 
(m, ri2, n^). Also longitudinal plane (77) containing jet axis - 
event plane, transverse plane (77 x) and components of particle 
momentum p in transverse plane are shown in the system of 
coordinate. 



Historically sphericity and thrust are considered as 
main traditional variables characterized the geometry of 
final hadronic state. There are several additional collective 
observables derived from the sphericity tensor. In particu- 
lar, p? n and p are, respectively, the squares of particle 
transverse momenta in and out of the event plane (Fig[T|). 
One of the versions for definitions of collective variables 
via eigenvalues of sphericity tensor S a ^ are shown in the 
Table □ 
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Table 1. Event shape observables 



are uncorrelated in different experiments, and calculate 
a weighted average as 



Observable 



Definition 



sphericity (5") 3(A 2 + A 3 )/2A E 



Pout 



alignment (A) 
planarity (P) 



A 2 /iV 
A 3 /iV 
A3/A1 
A 3 /A 2 



According to definitions J2| and J2]) S 1 and T -t 
0.5 for events with large number of the secondary particles 
distributed isotropically in phase space; S — > and T — > 1 
if two separated and narrow collimated particle groups are 
produced in interaction process which are directed to the 
opposite sides. Magnitudes of components p? and p out 
depend on a choice of a reference plane in a longitudi- 
nal direction, i.e. a plane contained jet axis (plane 77 at 
FigJT]) while absolute value of square of particle transverse 
momentum with respect to the jet axis p\ = p? + p^ ut 
does not depend on a choice of this plane (FigJTJ) . The 
detailed review of the experimental results obtained for 
choice as (longitudinal) scattering plane the plane formed 
by the two beams and the trigger particle is presented, for 
example, in [24]. Events with greater alignment are char- 
acterized by smaller value of parameter A, for strongly 
planar events is valid the relationship Pel. 

Symbol Y denotes below the generic event shape vari- 
able from the set under study. The moment of the 
distribution of an event shape parameter Y is defined by 



(Y n ) = 



dYY r - 



. 1 da 



n = l,2,... 



(5) 



Here Y miiX is the kinematically allowed upper limit of vari- 
able Y, crtot ~ total hadronic cross section. The present 
paper focuses on the first moments of distributions of col- 
lective variables. 



3 Experimental data sample 

Experimental results for event shape variables are con- 
sidered for all interaction types and for total available en- 
ergy range in the paper. One need to stress that although 
event shape of hadronic final state is studied for a long 
time there is no standard database for collective charac- 
teristics. The experimental data are presented in the HEP- 
DATA [25] and we made attempt to collect results for wide 
set of collective observables. Table [5] shows corresponding 
facilities and their experiments. Results from latest paper 
of the experiment are included for given collective param- 
eter at certain energy. 

To average data from different experiments at equal 
or, at least, (very) close energies weighted procedure is 
used. We assume that measurements of given quantity 



Y = 



Si 6* 



(6) 



where (Y)i is the average value of given observable re- 
ported by the i^ 1 experiment, £j is the weight of the 
i th experiment which is equal statistics / luminosity used 
for study of given observable in the experiment, and the 
sums run over the all experiments at certain energy. Ta- 
ble [3] shows energy values and corresponding experiments 
whose results were averaged for some given collective ob- 
servable. The notation ^/s is used for all interaction types 
but one need to stress that Q (absolute value of momen- 
tum transfer) and W (invariant mass of final hadronic 
state) imply as initial energy parameter for DIS experi- 
ments at HERA and for experiments with (anti)ncutrino 
beams respectively. 

It would be stressed that the common basis, in partic- 
ular, adequate initial energy scale needs for correct com- 
parison between various ways of producting multihadronic 
final states [74] . Identification of the correct initial energy 
scale, especially, for hadronic interactions, seems the non- 
trivial task. There are different methods for definition of 
energy scale in hadronic reactions (see, for example, [60, 
63, 66]) but universal scheme is absent so far. Detail stud- 
ies in the framework of the paper show that y/s can be 
used as an energy scale in hadronic collision together with 
some specific additional cuts for event / particle selec- 
tion (for example, reducing of influence of leading particle 
etc.). Comparison of average values of different collective 
observables will show reasonable agreement for various in- 
teractions at such choice of energy scale for hadronic in- 
teraction (see below). 

Also averaging has been made on Breit variable x, 
beam energy for some DIS experiments and on different 
measured kinematic ranges for some hadronic interactions 
in order to exclude the possible influence of additional 
kinematic cuts on the collective characteristic values. Er- 
rors of experimental points include available clear indi- 
cated systematic errors added in quadrature to statisti- 
cal ones. One need to emphasize that the most of data 
for main collective parameters (sphericity and thrust) are 
published results. On the other hand the some values of 
the above parameters as well as for other observables are 
estimations which were obtained on the basis of mean val- 
ues of eigenvalues of sphericity tensor or at experimental 
distributions of p- m , Pout- Obviously, such approach al- 
lows to get only rough estimations of true average values. 
But on the other hand the smooth experimental dependen- 
cies will be presented below and corresponding discussion 
show that estimations obtained by this method seem rea- 
sonable for any collective parameters under study and al- 
low to extend energy domain significantly. Also the values 
of main collective observables for special cases as well as, 
for example, sphericity and thrust values for quarkonium 
states [29, 30] or thrust values in the very narrow energy 
ranges at searching for top quark [36] are not shown at 
the global fits below. 
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t acuity 


Location 


Experiment / Detector 


Kets. 




1. e 4 


e~ annihilation 




SPEAR 


SLAC Stanford 


SLAC-LBL 


[13] 


PEP 




DELCO 


[26] 






HRS 


[27] 






Mark II 


[28] 


DORIS 


DESY Hamburg 


PLUTO 


[17, 29] 






Nal 


[30] 


PETRA 




CELLO 


[31] 






JADE 


[32 - 34] 






MARK-J 


[35, 36] 






PLUTO 


[17, 37] 






TASSO 


[38, 39] 


TRISTAN 


KEK Tsukuba 


AMY 


[21] 






TOPAZ 


[40] 


SLC 


SLAC Stanford 


Mark II 


[41] 






SLD 


[42] 


LEP 


CERN Geneva 


ALEPH 


[43] 






DELPHI 


[44] 






L3 


[45] 






OPAL 


[33, 46] 


2. interactions with lepton beams 


SPS 


CERN Geneva 


4m BC (BEBC) 


[47] 


PS 


FNAL Batavia 


4.6m BC 


[48] 


HERA 


DESY Hamburg 


HI 


[49 - 52] 






ZEUS 


[53] 


3. hadron-hadron collisions 


AGS 


BNL Upton 


2m BC 


[54] 


PS 


CERN Geneva 


0.8m and 2m BCs 


[54] 






2m BC 


[22] 






1.5m, 2m, and 4m BCs 


[55] 






4m BC 


[56] 


ISR 




CMOR 


[57] 






NA5 


[58] 






NA23 


[59] 






SFM 


[60] 


SppS 




UA2 


[61] 


LHC 




CMS 


[62] 


PS 


FNAL Batavia 


IHSC 


[63] 






LAFMS 


[64] 


Tevatron 




CDF 


[65] 


U-70 


IHEP Protvino 


2m BC JINR ("Ludmila") 


[23, 66] 






2m BC 


[67] 






4.5m BC ("Mirabelle") 


[68, 69] 


4. hadron-nuclear interactions 


PS 


ITEP Moscow 


lm and 2m BCs 


[70, 71] 


U-70 


IHEP Protvino 


2m BC JINR 


[72] 


SPS 


CERN Geneva 


BEBC 


[47] 


PS 


FNAL Batavia 


E260 


[73] 
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Table 3. Energy and corresponding experiments for averaging 



y/s, GeV 


Experiments 


1. 


e + e~ annihilation 


9.4 


Nal, PLUTO 


12.0; 13.0; 17.0; 27.3; 30.6 


MARK-J, PLUTO, TASSO 


14.0 


JADE, TASSO 


22.0 


JADE, MARK-J, PLUTO, TASSO 


29.0 


DELCO, HRS, Mark II 


35.8 


CELLO, JADE, MARK-J, TASSO 


38.3 


JADE, MARK-J 


43.6 


CELLO, JADE, TASSO 


M z 


ALEPH, DELPHI, L3, Mark II, OPAL, SLD 


133.1 - 206.1 


ALEPH, DELPHI, L3, OPAL 


2. interactions with lepton beams 


3.2; 3.9; 5.0; 6.5; 10.6 


BEBC, FNAL BC 


9.3; 13.3; 18.3 


HI, ZEUS 


11.0; 15.2; 17.7; 22.6; 29.7 




36.8; 42.6; 58.9; 82.3; 116.3 





As known definitions of collective observable may vary 
from one experiment to other and depend on interaction 
type and energy domain under study (for example, there 
are several definitions for thrust analogies). But in the 
framework of this paper we have tried (if it was possibly) 
either to estimate the average values of collective param- 
eter or recalculate they in according with the definitions 
shown in Table [T] Moreover, for thrust we chose in original 
papers the definition of the variable which is most close 
to the (J2J) both in mathematical and physical meaning for 
all interactions in high energy domain. 



4 Energy dependence 

On Fig|3] dependencies of average values of collective 
observable sphericity (a) and thrust (b) on collision en- 
ergy are shown. Experimental results and the estimations 
for collective characteristics under study are taken from a 
database presented in Table [2 As seen from Fig J2] general 
tendency for interactions of various types is observed both 
for sphericity (a) and for thrust (b) at a qualitative level, 
moreover the discrepancy of experimental points reduces 
at increasing of y/s. 

In the case of hadronic reaction influence of in the 
leading particles results in decreasing (increasing) of av- 
erage value of S (T) in comparison with e + e~ and IN 
collisions. Selection of events with high multiplicity, spe- 
cific type of interaction (non-diffractive) or using of mod- 
ified sphericity tensor allows to reduce the influence of 
the leading particles on jet characteristics in hadronic re- 
actions. Additional cuts indicated above result in better 
agreement of values of (S) and (T) in various interactions 
in the range of intermediate energies. Experimental re- 
sults for 7r~Nc collisions at 6.2 GeV/c [70, 71] demon- 
strate the influence of nuclear matter on jet observables. 



Exception of interactions with single intranuclear nucleon 
and, probably, small, in comparison with radius of a tar- 
get nucleus, length of forming of final hadronic state [75] 
leads to essential influence of multiple rescattering on jet 
characteristics which results in wider jets and significantly 
larger (smaller) average values of sphericity (thrust) in 
comparison with results of other experiments at close en- 
ergies. The account of interactions with a single intranu- 
clear nucleon in case of discussed 7r~Ne interactions leads 
to (S) = 0.269±0.005 and (T) = 0.814±0.002, that would 
be agreed with the general tendency much more reason- 
ably. Reducing of influence of rescattering with increasing 
of beam momentum leads to improvement of agreement 
of (S) even for multinucleon hadron-nuclear interactions 
with the general tendency (Fig|3k), that is visible on an 
example of ir~C collisions [72]. Thus, integral (on num- 
ber of participating nucleons) results confirm a conclusion 
[71, 72], that the mechanism of production of hadron jets 
at nuclei is more complex, than that in processes with 
hadrons and / or leptons, and growth of average values of 
sphericity for hA interactions at intermediate energies can 
be explained by multiple rescattering of secondary parti- 
cles at nucleons of a target nucleus. Point for 7r~Ne at 
initial momentum 6.2 GeV/c discussed above is excluded 
below from the fitted samples for sphericity and thrust. At 
present there are no experimental results in TeV energy 
domain for traditional collective observables which were 
described above. The experimental results for transverse 
thrust (Tj_) from [62, 65] are used as estimations of thrust 
defined by ([3]) at TeV scale energies. This assumption is 
based on the relation T ~ T±_ which valids for events with 
two back-to-back jets and for near-to-planar 3-jet events 
[76, 77, 78]. As seen from FigJ2]the CDF experimental re- 
sult [65] differs noticeably from general tendency for thrust 
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■ 7t-(C 2 F 5 CI 3 ) ■ h + p X K + p O ep _ 

□ 7t"Ne A 7t"p ▼ PP • (v+v)N 

■ nC O nn >k K'p e + e" 

▲ PP 



Fig. 2. Average values of the sphericity (a) and thrust (b) depend on initial energy. Experimental results and estimations 
from database presented in Tabled Total errors are shown. See text for explanation of the curves. 
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Table 4. Results of fit of dependence (S)(s) by function (J7J 



Data sample Fit parameters 







a 2 


a.3 


X 2 /ndf 


hh and hA 


0.047 ± 0.004 


0.442 ± 0.005 


-0.210 ±0.006 


97.5 


IN 


-357 ± 303 


358 ± 303 


(-1.5 ± 1.2) x 10" 4 


5.22 


e+e" (y^ > 5 GeV) 


0.019 ±0.004 


0.94 ± 0.02 


-0.311 ±0.009 


2.95 


all data 


0.021 ±0.002 


0.610 ± 0.004 


-0.244 ± 0.004 


70.0 



observabl^E i.e. the assumption is less correct for CDF 
than that for CMS. Possibly there are several reasons for 
such behavior of experimental point at y/s = 1.96 TeV. As 
known contribution of soft particles at large angles with 
respect to the event plane can leads to T > T± [77]. In- 
deed, CDF data show slightly intensive energy flow out of 
the hard-scattering plane [65]. Thus discrepancy of CDF 
point from general trend can be caused, in particular, by 
more significant contribution of soft particles at large an- 
gles. Additional reason may be kinematical cuts used in 
CDF analysis. In any case for verification of this hypoth- 
esis additional detail study is essential for event geometry 
with taking into account soft particles, especially, exami- 
nation of relation p m ~ Pout f° r such particles. Therefore 
the CDF experimental point is excluded below from fitted 
samples. 

There are smooth dependencies of experimental (S) on 
y/s for various interactions without qualitative changes at 
transition from intermediate energy domain to high ener- 
gies (Fig[2k) • Taking into account general view of spheric- 
ity dependence on initial energy and results for e + e~ [27, 
39], in the present paper the (S)(s) is approximated by 
power function 

(Y) = ai + a 2 X a -\ (7) 

where Y is certain collective variable and X denotes the 
appropriate argument of a function, Y = S and X = s in 
the case under considered. Fit results for various samples 
of experimental data for sphericity observable by func- 
tion Q are presented in Tableland shown on FigJ5^ for 
hh ± hA interactions by a solid line, for IN - by dashed 
one, for e + e~ annihilation - by dotted line and for the 
joint sample by thin solid one. The disorder of experimen- 
tal points, especially in the range of \fs < 5—10 GeV, does 
not allow to obtain statistically acceptable fit quality. In- 
fluence of discrepancy of experimental data from indicated 
above energy domain on the approximation quality is con- 
firmed, for example, by results for e + e~ annihilation: for 
y/s > 5 GeV reasonable fit quality is observed (Table 2]), 
moreover value of 03 parameter agrees within errors with 
result from [27], on the other hand, fitting of all available 
energy range results in x 2 /ndf = 19.1. However one can 
see from Fig[2l that the behavior of (S) depending on yfs 
is described by function ([7]) at a qualitative level both for 

1 One need to emphasize that the result based on parton- 
level NLO±NLL distribution agree much better with general 
trend. 



separate interactions, and for the joint sample for all avail- 
able range of initial energies. Additional research for the 
joint sample with using of more complex fit function, in 
particular, the sum of power functions, has shown, that 
experimental data are better described by function ([7]). 

Unlike sphericity, collective observable the thrust pos- 
sesses infrared stability that allows to use QCD for study 
of behavior various dependencies and distributions on this 
collective parameter. In the framework of pQCD the en- 
ergy dependence of the first moment of thrust distribution 
can be described by the following function 

( T )pQ cd ( s ) = 1-£a4(s), <*s{s) = (bot)- 1 . (8) 

i=l 

Here t = \ns/A 2 , A is the QCD parameter, bo is referred 
to as one- loop /3- function coefficient [2]. Parameters Ai 
are calculated up to NNLO in [79] and for ([8]) they are 
Ax = 0.335, A 2 = 1.033, A3 = 2.76 ± 0.09 for num- 
ber of (light, i.e. active) quark flavors Nf — 5. FigJ^b 
shows the theoretical predictions for LO ("1"), NLO ("2") 
and NNLO ("3") pQCD based on the © by thin dotted 
curves. As seen behavior of a theoretical curve (T) pC ^ CD (s) 
depends on the order of pQCD expansion, especially, in 
the intermediate energy domain. The agreement of theo- 
retical prediction (JSj) with experimental data improves for 
higher order of pQCD expansion. The LO pQCD overpre- 
dicts the thrust value even at LHC energy. It is important 
to note very good agreement of experimental results for 
pp collisions at y/s = 7 TeV [62] and a theoretical curves 
(© for NLO and NNLO. Such quantitative agreement is 
observed for the first time, especially, for hadronic reac- 
tion and allows to assume, that at yfs ^5 — 10 TeV non- 
perturbative corrections become negligible. As known the 
contribution of these corrections increases at decreasing 
of initial energy and for approximation of (T)(s) in range 
of yfs as wide as possible it is necessary to use the phe- 
nomenological models which take into account the non- 
perturbative corrections. Thus, in the present work for fit 
of dependencies (T)(s) in various interactions the follow- 
ing general function is used: 

(T)(s) = l-[<F) pQCD + (F) n -P], (9) 

where components (F) p Q CD and (F) n ~ p are model depen- 
dent. In the framework of present paper dispersive model 
(DM) [80, 81] and single dressed gluon model (SDGM) 
[82 - 84] have been used for approximation of dependence 
(T)(s). It is necessary to note, that now the most part 
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of calculations are executed for electron-positron annihi- 
lation and comparison of models with experimental val- 
ues of collective characteristics is carried out, as usual, 
for the this type of interactions because precision of mea- 
surements and calculations for shape variables of final 
hadronic state is the best for e + e _ namely, in additional 
to significant sample of experimental results is available in 
the high energy domain where applicability of QCD and 
corresponding phenomenological models has more rigor- 
ous substantiation. Some predictions on the basis of DM 
have been obtained for ep interactions also [49, 51, 52]. For 
hadron-hadron and hadron-nuclear reactions the number 
of phenomenological predictions is much less, that is con- 
nected both with the raised difficulty of calculations and 
with absence of experimental data at high energies until 
recently. Therefore the functional forms obtained for e + e~ 
initially are used below for fitting of various data samples. 

In DM perturbative part is defined in the framework 
of NLO pQCD and terms for are the following: 

2 



(F) 



n-p _ A// 4Cf f flj 
DM 



K 



(10) 



a (/i/) - a s {s) - b In — + — — + 2 a s (s) 
\ [ij 2irb J 

V. 18 6 ) 9 11 3' 



Here Ai is correcting so called Milan factor [85, 86]. In 
this work M. and Ai, (i — 1,2) are free fit parameters. 
Non-perturbative parameter 



p + l 



k p a s (k)dk 



is introduced for the account of a part of expression for 
as(s) which diverges below fii - infrared matching scale 
[80] . There is no unambiguous scheme for definition of fii 
value. The empirical choice of fij and Nf values is re- 
stricted by the requirement of validity of following condi- 
tions: A <C Mi < V s min [80] and ViV/ : m q < ^/s mm [2], 
where y/s m i n is the low boundary of fitted energy range, 
m q denotes masses of the quarks, which can be considered 
as lights for a certain ^/s mm . One need to emphasize that 
in the framework of present paper, in distinction from, for 
example [51], for IN collisions dependencies of parameters 
Ai, (i — 1,2) on some kinematic quantities (x, Q) have 
not been taken into account because, as indicated above, 
the averaged on these parameters and various experiments 
values of thrust are used. Fit results for various samples of 
experimental data for thrust observable by function @ in 
the framework of DM are presented in Tableland shown 
on FigJ2h> for e + e~ annihilation by dotted curve, for IN - 
by dashed one, for hh+hA interactions - by solid line, and 
for the joint sample - by thin solid one. As usual the bet- 
ter fit quality is obtained at the fixed value of M which is 
defined by limits of the allowable interval of Milan factor 



values [85]. As seen from Table [S] and FigJ2b DM allows 
to describe all available results for e + e~ at reasonable fit 
quality. The disorder of experimental points leads to that 
statistically acceptable values of x 2 /ndf are observed at 
^/s mm = 5 and 10 GeV only for IN and hadronic reactions 
respectively. For these interactions the fit functions are 
shown Fig|2jD at ^/s mm indicated above namely. However 
function © in the framework of DM allows to describe 
dependence (T)(s) at a qualitative level both for some in- 
teractions and for the joint data sample in all range of 
available initial energies (FigJ^b) . Behavior of DM-curves 
agrees very well with NLO and NNLO pQCD predictions 
in the LHC energy domain. Thus, for the first time the 
reasonable phenomenological description is obtained for 
experimental results for thrust variable in all available en- 
ergy range both for e + e~ and for the joint data sample for 
various interactions. Non-perturbative parameter a p (jii) 
values obtained for e + e and IN at y/s m [ n = 10 GeV 
agree with previous results [52, 87] within errors. Clear 
dependence of a p (/j,i) on ^/s m j n is absent within errors 
for these interactions. The similar situation is observed 
for hadronic reactions. But one can see growth of a p (/i/) 
value at decreasing of ^/s m j n for joint data sample, es- 
pecially visible at transition in strong non-perturbative 
region from y / s m ^ n = 5 GeV down to ^/s mm = 2 GeV 
(Table O). There is no clear dependence of a p (ni) value 
on type of interaction at fixed low boundary for fitted en- 
ergy domain. 

For e + e~ results and for joint data sample approxi- 
mations in the framework of SDG-model have been in- 
vestigated too. In the case of the first moment of thrust 
distribution, (T) , perturbative and non-perturbative com- 
ponents in SDGM for © are the following [87]: 

<^>S§GM=I>«W^ ( f )sDGM = ^- (11) 



Here 



ais 



a s {s) 



1 - ^b a s {s) 



and v is a free fit parameter. Coefficients Bi are defined 
via so called log moments di of the characteristic thrust 
function, calculated in [83] , by following way 

8i = do, B 2 = pod u B z=( y -^- + $ + & d i> 
B 5 = U + - 2^ 2 d 2 ) /3 4 + (Hftrf, _ fr^A ft 

+ 3/3^2/30 + ^2+/^!, 

B 6 = f-y^da + ^di + (kj /3 5 + (^d 4 - -tt 2 ^ x 

/3i/3 3 + (6d 3 - **di)fafi + (y/Sids + \kd^j A)+ 
7 

Mi + -^ifodi- 
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Table 5. Results for fit of dependence (T)(s) by function (|9} in the framework of DM 



Fit 


Data sample 




parameter 


e+e" IN hh and hA 


all data 


= 2 GeV (m = 1 GeV, N f = 3) 



Ai 3.8 ±0.4 2.1 ±0.3 4.08 ±0.15 2.6 ± 0.3 

Ai -63 ± 26 48 ±17 -11.1 ±0.8 16 ± 6 

M (fixed) 1.795 1.448 1.100 1.100 

aoQn) 0.77 ±0.04 0.57 ±0.03 0.628 ± 0.012 0.75 ±0.05 

as(Mz) 0.1192 ±0.0007 0.1209 ± 0.0012 0.1180 ± 0.0011 0.1218 ± 0.0007 

X 2 / ndf 3.04 8.16 20.1 15.6 



ysT^ = 5 GeV Qxi = 2 GeV, Nf = 4) 



Ai 3.17 ±0.13 8.02 ±0.06 4.2 ± 1.2 2.76 ±0.15 
A 2 -32 ±4 -300 ±7 6±3 2.7 ±1.3 
M (fixed) 1.520 1.852 1.188 1.188 
So (A*/) 0.532 ±0.018 0.638 ± 0.008 0.39 ± 0.02 0.517 ± 0.008 
a s {M z ) 0.1203 ±0.0018 0.1222 ± 0.0003 0.116 ±0.004 0.1206 ± 0.0009 
X 2 /ndf 3^63 L74 203 7.38 



^/s^ = 10 GeV (jmi = 2 GeV, 7V/ = 5) 



Ai 4.4 ±1.3 -0.6 ±0.3 2.7 ±1.1 3.7 ±1.2 

A 2 -97 ± 66 212 ±21 8±2 -58 ±34 

M 1.291 (fixed) 1.291 (fixed) 1.920 (fixed) 1.8 ± 0.5 

ato{fJ,i) 0.59 ±0.18 0.32 ±0.17 0.42 ±0.10 0.49 ± 0.04 

a s {M z ) 0.115 ±0.011 0.115 ±0.017 0.117 ±0.016 0.121 ±0.011 

X 2 /ndf 3.28 0.43 0.69 3.19 



Table 6. Results for fit of dependence (T)(s) by function in the framework of SDGM 



Data 


Fit parameter 




sample 


a s (M z ) v 


X 2 /ndf 


e + e~ 


0.1235 ±0.0008 -0.28 ±0.09 


6.98 


all data 


0.1197 ±0.0005 -0.142 ±0.049 


14.0 



Here /3j — (4iry +1 bj, j = — 4, bj are coefficients of 
/3-functions in j-loop approach [2] with that /?4 = be- 
cause of 64 is unknown [87]. Scheme described above for 
DM is used in this case too for choice of Nf depending on 
y /s m i n value. Detail study shown, that at the fixed val- 
ues of di from [83] SDGM allows to obtain the description 
of experimental data at a qualitative level for both con- 
sidered samples only at ^/s mm = 10 GeV (Table [BJ with 
that presence of non-perturbative term weakly influences 
at fit quality. Expansion of fitted range down to the lower 
energies leads to significant discrepancy of fit curve with 
experimental data. 

Fig (3] shows values of as(Mz) obtained by fit of de- 
pendence (T)(s) in the framework of DM and SDGM at 
various ^/s mm . Experimental values of as(Mz) derived 
a t v^min = ^ and 5 GeV for some interactions as well 
as for joint data sample agree reasonably with world av- 
erage [2, 88]. Precision for estimations of strong coupling 
constant considerably deteriorates at ^/s mm = 10 GeV 
in comparison with smaller values of the low boundary of 
fitted energy range in the framework of DM (Fig|3fc). At 



this value of ^/s mm SDGM allows to improve significantly 
the precision of as{Mz) estimation. The value of as(Mz) 
extracted in the framework of SDGM for the joint sample 
is in good agreement with world average value Og a (Mz). 
On the other hand there is some excess of experimental 
value of strong coupling constant for e + e~ sample with 
respect to the dg & {Mz)- However as(Mz) experimental 
value shown in Table |5] for e + e~ coincides within errors 
with the results obtained, for example, at investigation of 
radiative decays of heavy quarkonium (T) and shape of 
final hadronic state in ALEPH experiment [2]. Thus, in- 
vestigation of energy dependence of average thrust values 
for various interactions allows to extract estimations of 
cts(Mz) which agree reasonably both with world average 
and with results of other methods. 

Figure [4] shows energy dependence for (pf n ) (a), (p ou ^) 
(b), (A) (c) and (P) (d). Experimental results and the 
estimations for collective observables are taken from a 
database presented in Table [2 Dependencies for all collec- 
tive parameters under study show smooth behavior, with- 
out qualitative changes at transition from intermediate 
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e + e" 



IN 

hh+hA 
all data 



e + e" 



IN 

hh+hA 
all data 



e + e" 



IN 

hh+hA 
all data 



\ s . = 2 GeV 

\ mm 



\ s . = 5 GeV 

\ mm 




\s min = 10 GeV 

\ mm 



□ 



0.1 



0.11 0.12 

a s (M ) 



0.13 



Fig. 3. Values of as(Mz) obtained by fitting of various samples of experimental data at ^s mm = 2 (a), 5 (b) and 10 GeV (c). 
Values extracted in the framework of DM are shown by closed symbols, for SDGM - by opened ones. Dashed line and shared 
band are shown the world average value with errors as & {M z ) = 0.1184 ± 0.0007 [2, 88]. 



energy domain to high energies (FigU). Unlike the situ- 
ation for one of the main collective parameters, T, the 
amount of phenomenological models describing behavior 
of observables depending on initial energy considered on 
FigQJ is essential less or they are absent at all. There- 
fore in the framework of present paper the wide set of 
functions has been studied for empirical approximations 
of experimental results for total available energy range. 

Experimental results for hadron- nuclear 7r~Ne interac- 
tions agree both with results for hadron-hadron reactions 
at similar energies and with general trends reasonably in 



case of a component of a transverse momentum in "prin- 
cipal axes" frame (FigHk and FigQJ)). On the other hand 
the multinucleon 7r~Nc collisions show essentially larger 
values of parameters (A) (Fig|4];) and (P) (Figj4jl), i.e. 
hadron-nuclear collisions at y/s ~ 6 GeV are character- 
ized by weaker alignment of hadrons with respect to the 
jet axis and smaller planarity of events in comparison with 
hadronic interactions, that can be explained, possibly, by 
rescattering of hadronic states in nuclear matter discussed 
above. Therefore the experimental results for 7r _ Ne are ex- 



12 



V.A. Okorokov: Event shapes for hadronic final state: experimental review 




Fig. 4. Average values of the parameters describing the event shape in the "principal axes" coordinate system depend on initial 
energy: a - (p m ), b - (p ou ^), c - (A), d - (P). Values of (p? n ) and (Pq U ^) are shown in (GeV/c) 2 . Experimental results and 
estimations are from database presented in Table [2] Total errors are shown. See text for explanation of the curves. 
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Table 7. Fit results for energy dependence of collective observables under considered 



Data 


Function 




Fit parameter 






sample 




ai 






X 2 /ndf 


(Pt) 


hh and hA 


G3 


-36 ± 28 


36 ±28 


(1.8 ± 1.4) x 10 -3 


51 


ep 


-II- 


-0.08 ± 0.03 


0.06 ± 0.02 


0.7 ±0.3 


0.53 


e + e~ 


-II- 


0.033 ±0.025 


(1.6 ±0.8) x 10~ 3 


2.5 ±0.3 


2.24 


(Pint) 


hh and hA 


G2J 


0.007 ± 0.003 


(2.1 ±0.8) x 10~ 3 


1.5 ±0.3 


31 


ep 


-II- 


-0.119 ±0.009 


0.099 ± 0.007 


0.26 ± 0.02 


0.19 


(hh, hA)+ep 


-II- 


(1.08 ±0.02) x 10 -2 


(1.6 ±0.04) x 10" 3 


1.0 (fixed) 


12 


e + e~ 


-II- 


0.024 ± 0.007 


(4± 2) x 10" 4 


2.4 ± 0.3 


0.79 


(A) 


hh 


(7) 


-0.19 ±0.07 


0.57 ±0.06 


-0.12 ±0.03 


195 


e + e~ 


-II- 


0.027 ±0.004 


0.7 ±0.2 


-0.42 ±0.06 


2.29 


all data 


-II- 


0.017 ±0.002 


0.483 ± 0.008 


-0.323 ± 0.008 


25.3 


<P> 


ep 


G3J 


0.30 ±0.02 


(2.6 ± 1.4) x 10~ 4 


(-2.8 ± 1.7) x 10" 7 


0.16 


e + e~ 




-3.24 ±0.07 


3.85 ±0.07 


-0.012 ±0.001 


0.88 



eluded from the fit samples in the case of (A) (Fig0t) and 
(P) (FigHH) variables below. 

Comparison of FigHk with FigH)D demonstrates that 
value of " in" -component of transverse momentum is larger 
than corresponding "out" -component value at all consid- 
ered initial energy values yfs > 4 GeV and for all inter- 
action types. In the energy domain accessible for compar- 
ison, values of (pf n ) (FigH^) and (p^) (Fig0b) show 
faster growth at increasing of y/s in case of e + e~ annihi- 
lation in comparison with interactions of other types. It is 
agreed with results of the comparative analysis executed 
in [60] for mean square of transverse momentum, (p 
at 10 < y/s < 35 GeV. In the framework of present pa- 
per energy dependence for first moments of distributions 
of "in"- and "out" -component of transverse momentum is 
fitted by function 

(p 2 J )(s)=a 1 +a 2 [lns} a *, j = in, out. (12) 

Here s is in GeV 2 . Results of fit of components of trans- 
verse momentum by function (|12[) are presented in Table 
[7J for various experimental data samples and shown on 
Fig|4^,b for hh and hA collisions by a solid line, for ep 
- by dashed line and for e + e~ annihilation - by dotted 
one. The behavior of (pf n )(s) and (p ou t)( s ) depends on 
type of interaction, especially visible for the first observ- 
able. Disorder of experimental points and small errors do 
not allow to get reasonable fit quality for hadronic interac- 
tions, however function (|12l) agrees with the general ten- 
dency of experimental results at a qualitative level. Usu- 
ally experimental data and estimations for e + e~ and ep 
are described by the offered function (fT2"j) with statistically 
acceptable quality both for (p m ) (FigHk), and for (Po Ut ) 
(FigHJ}). As it has been stressed above, the estimations 



derived on the basis of published e + e~ experimental data 
are used for, in particular, (p m ) and (j> ou t) a t yfs > Mz- 
Direct experimental measurements are available for these 
observables up to y/s — 55.2 GeV in e + e~ annihilation 
[21]. Therefore results for e + e _ have been studied in en- 
ergy range yfs < 60 GeV separately, in which (pfL) and 
(Pq U ^) obtained from experiment data directly. For pa- 
rameter (pj n ) better fit quality turns out for power func- 
tion ([7]) at Y = p? in case under considered. Fit func- 
tion for energy range indicated above is shown on Fig|3^ 
by thin dotted line, values of fit parameters are follow- 
ing: ai = 0.124 ± 0.005, a 2 = (9.3 ± 0.6) x 10" 5 , a 2 = 
1.0 (fixed), x 2 / n df = 1-62. Thus, measurements only at 
y/s < 60 GeV assume significantly faster growth of (p ? ) , 
than that observed at the account of estimations at LEP 
energies. For (j>g Ut ) the published results and estimations 
at collision energies larger than Mz show very smooth 
behavior without any features. Values of fit parameter 
for (|12l) are coincide within errors for fitting of all en- 
ergy range and for -^/s < 60 GeV domain only. Fitting by 
function (IT2l has been made for joint sample of results 
for hadronic and ep interactions because of close values 
of (j5Q Ut ) for these collisions. The values of fit parameters 
presented in Table [7J close within errors with results for 
hh ± hA data sample. 

Alignment of secondary particles with respect to the 
jet axis amplifies with increasing of collision energy both in 
hadronic and in e + e~ interactions (FigHJ:). In the frame- 
work of present paper dependence (A) (s) is fitted by func- 
tion ([TJ, where Y = A in case under study. Fit results for 
various samples of experimental data are presented Table 
[TJand shown on FigHJ: by solid line for hh interactions, for 
e + e _ annihilation - by dotted one and for sample of all 
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available data - by thin solid curve. As seen from FigHJ; 
the behavior of collective observable (A) depending on 
initial energy is described by function (|7J| for all sample 
under study at a qualitative level. 

Absence of experimental data in the range of high / 
low energies for (hh, ep) / e + e~ interactions respectively 
allows to propose only qualitative hypothesis concerning 
of behavior of planarity parameter depending on initial 
energy. Observable growth of (P) up to y/s ~ 30 GeV 
evidences on broadening jets in a transverse plane, that, 
probably, is dominated by the contribution of soft particles 
and more isotropic distributions of particle momenta in a 
transverse plane. At increasing of initial energy events be- 
come more and more planar (FigHJi). In the present work 
dependence of planarity parameter on collision energy for 
ep interactions is approximated by a polynomial function 



JY 



<p> = £ 



dfcS 



(13) 



k=0 



of 2-d order and by function (0 - for e + e~ annihilation. 
Results for fitting of these experimental data samples are 
presented in Table [7] and shown on Fig (Hi by dashed curve 
for ep interactions and by dotted line for e + e~ annihila- 
tion. As fitting function of sample of experimental data for 
hh interactions had been considered power function (|7|), 
predicting smooth growth of planarity in collisions of the 
specified type at high energies, and polynomial function 
(|13p of various orders (in the range of N = 2 — 8). Better 
quality of approximation among the studied functions is 
observed for function (fT3")) at N = 5 (x 2 /ndf = 11.9) for 
which following values of parameters have been obtained: 
at = 0.192 ± 0.002, a 2 = (1.84 ± 0.07) x 10" 3 , a 3 = 
(-9.1 ± 0.5) x 10- 6 , a 4 = (1.71 ± 0.11) x 10" 8 , a 5 = 
(-1.05 ± 0.08) x 10 -u . This fit by function (13) for hh 
interaction sample shown on FigHJl by solid line, demon- 
strate fast decreasing of (P) at y/s > 30 GeV. The reason 
of such behavior, as well as in case of ep collisions, pos- 
sibly, is absence of experimental data at y/s > 30 GeV. 
Thus, additional experimental data are necessary for the 
quantitative analysis and more unambiguous conclusion 
concerning the behavior of (P)(s) in hadronic interactions 
at high energies. During investigation of (P) for all domain 
y/s > 2 GeV different versions of additional term in (J7]) 
have been examined, in particular, Landau, normal and 
log normal functions. These researches have shown, that 
joint sample of experimental results for hh, ep and e + e _ 
collisions supposes fit by the function 



(P) = ai +a 2 s° 3 + — exp 
s 



1 / In s — a§ 

2 ^ a~ 6 



(14) 



with the best x 2 / n df = 5.36 among studied functions 
and with the following values of fit parameters: a\ — 
0.11 ± 0.02, a 2 = 2.3 ± 1.3, a 3 = -1.8 ± 0.9, a 4 = 
-1043 ± 587, a 5 = 11.2 ± 0.9, a 6 = 2.3 ± 0.2. As seen 
from FigHH the suggested approximation function shown 
by thin solid line agrees qualitatively with experimental 



results at intermediate energies and corresponds to ex- 
pected amplification of event planarity in high energy do- 
main. 

On the basis of the above experimental results, esti- 
mations and the suggested approximations it is possible 
to draw an intermediate conclusion, that all traditional 
collective variables show smooth dependence on energy 
without any features at y/s ^3 — 5 GeV. In Table [5] es- 
timations are shown for parameters of event geometry in 
" principal axes" coordinate system in various interactions 
for the present and possible future accelerating facilities, 
calculated in the framework of the paper on the basis of 
analytical functions (ffjl. ([9]) - ([TT ]) . (|TS |) - (jl"4)) and Table 
2]- EJ Additional experimental results are essential for im- 
proving of precision of estimations for (p m ) and (povrfc)- 
But it would be emphasize, that the estimation of abso- 
lute value of transverse particle momentum with respect 
to the jet axis y/ (p\) = 0.50 ± 0.17 GeV/c derived in the 
framework of present paper for pp collisions at y/s = 200 
GeV coincides within errors with experimental result at 
the same energy [89]. Study the high and ultra high en- 
ergy domain show that the asymptotic values within errors 
are predicted for most collective parameters under study 
at y/s ~ 5 - 7 TeV already. 



5 Multiplicity dependence 

Unlike energy dependence, the investigations of geome- 
try of a hadronic final state at the fixed secondary par- 
ticle multiplicities are less extensive and have been made 
for energy range y/s < 60 GeV. Hadron-hadron interac- 
tions were studied in the exclusive approach in the some 
cases for intermediate energies. Results for collective ob- 
servables were obtained only for jets of all identified sec- 
ondary particles in some experiments (see, for example, 
[22]). Moreover using of total multiplicity (iVtot) allows to 
increase amount of experimental points and hence to in- 
vestigate behavior of dependence of jet characteristics on 
multiplicity in more detail. Therefore the event shape of 
a hadronic final state is studied depending on multiplicity 
of identified particles (N), which is equal N = A^ t for 
the exclusive approach and N — iV cn in other cases. 

Average values of main collective variables are pre- 
sented at FigJS] depend on multiplicity of the identified 
secondary particles in event. Experimental data are taken 
from [22, 23, 48, 56, 63, 67, 69 - 71, 90]. At (W) > 4.0 GeV 
right points correspond to integrated values (S) for a mul- 
tiplicity range with low boundary indicated for (W) = 
4.0 GeV namely ((5^)- We used the results for modified 
sphericity tensor with reduced momenta of final state par- 
ticles [16] and sphericity definition from original papers 
for K~p [22, 69]. Fit results of present work and corre- 
sponding curves at FigJS] are described below. The inten- 
sification of spherical event shape is observed with final 
state multiplicity increasing for all types of interactions 
under study at intermediate energy domain y/s < 60 GeV 
(FigJS). Rescattering of hadronic states in nuclear matter 
results in to the some increasing of (S) for 7r _ Ne in com- 
parison with elementary reactions. Therefore one can see 
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Table 8. Estimations for event shape observables* 











Observable 






TeV 


(S) 


<T) 


(pi), (GeV/c) 2 {p 2 „ t ), (GeV/c) 2 


(A) 


<P> 








e + e 


annihilation 






0.5 


0.050 ±0.002 


0.955 ± 0.005 


0.9 ±0.8 


0.20 ±0.19 


0.026 ± 0.002 


0.11 ±0.02 


0.8 


0.044 ±0.002 


0.959 ± 0.005 


1.1 ±0.9 


0.23 ±0.23 


0.024 ± 0.002 


-II- 


1.0 


0.042 ± 0.002 


0.960 ± 0.005 


1.2 ± 1.0 


0.25 ±0.25 


0.023 ± 0.002 


-II- 


3.0 


0.033 ± 0.002 


0.966 ± 0.004 


1.7 ± 1.5 


0.35 ±0.37 


0.020 ± 0.002 


-II- 


5.0 


0.031 ±0.002 


0.968 ± 0.004 


2.0 ± 1.8 


0.40 ± 0.44 


0.019 ±0.002 


-II- 


ep interactions 


0.8 






0.25 ±0.25 


0.08 ±0.02 






1.6 


0.038 ± 0.002 


0.963 ± 0.005 


0.28 ±0.27 


0.08 ±0.02 


0.021 ±0.002 


0.11 ±0.02 








PP 


collisions 






0.2 


0.067 ±0.003 


0.946 ± 0.007 


0.17±0.16 


0.08 ±0.05 


0.033 ± 0.002 


0.14 ±0.03 


0.5 






0.18 ±0.18 


0.10 ±0.07 






7.0 


0.029 ± 0.002 


0.969 ± 0.004 


0.21 ±0.20 


0.16 ±0.13 


0.019 ±0.002 


0.11 ± 0.02 


14 


0.027 ± 0.002 


0.971 ± 0.004 


0.21 ±0.21 


0.17 ±0.15 


0.018 ±0.002 


-II- 


28 


0.025 ±0.002 


0.973 ± 0.003 


0.22 ±0.21 


0.19 ±0.17 


-II- 


-II- 


42 


0.024 ±0.002 


0.974 ± 0.003 


-II- 


0.20 ±0.18 


-II- 


-II- 


100 


0.023 ± 0.002 


0.976 ± 0.003 


0.22 ±0.22 


0.23 ±0.21 


-II- 


-II- 


200 


-II- 


0.977 ± 0.003 


0.23 ±0.22 


0.25 ±0.23 


-II- 


-II- 



*Fit results for joint data samples are used to estimate the (S), (T), (A) and (P) for a certain type of interaction, 
thus the some estimation is indicated once if various interactions are studied at the same energy. 



the good agreement of results obtained for 7r~Ne and for 
[v+v)y> in despite of some difference in the initial energies. 
Similar accordance with general tendency is observed for 
7T~p, K + p, pp and (v+v)p interactions at larger initial en- 
ergies too. In [15] analytical dependence (S)(N) oc N~ a - 5 
was derived for isotropic distribution of secondary parti- 
cles in phase space for two ultimate cases: non-relativistic 
and extremely relativistic ones which are shown at FigJS^i 
by dotted thin curves " 1" and " 2" respectively. Results for 
7r~Ne and for (y ± v)\> agree well with phenomenological 
dependence for ultra-relativistic case, that confirms the 
difficulty of separation of jet events at ^fs ^3 — 5 GeV on 
the basis of one event shape observable © for interactions 
of any types. Thus, one can suggest, that the universal 
estimation of the low energy boundary for experimental 
appearance of event jet structure in multiparticle produc- 
tion processes is ^fs~ c ~ 3 GeV. At increasing of initial en- 
ergy one can see that estimations extracted for isotropic 
one-particle distribution are essentially larger than exper- 
imental values (S) , that corresponds to an amplifying for 
appearance of event jet structure. Functional behavior of 
dependence (T)(7V) in case of e + e~ annihilation at sig- 
nificantly larger initial energy qualitatively distinguishes 
from behavior of (T)(N) for hadron-hadron collisions at 
y/s < 12 GeV but (T) decreases with multiplicity increas- 
ing for all interaction types under study (Fig{5jj) . One 
need to emphasize that universality of geometry of events 
can be considered as evidence for similarity of dynamics 
of soft hadronic jet production in various interactions in 
wide initial energy range. Therefore dependencies (S)(N) 
and (T) ( N) indicate at possible experimental appearance 



of color degrees of freedom in hadronic jet production at 
y/s - 3 - 5 GeV. 

Figure [6] shows the multiplicity dependence for average 
values of a square of in-plane transverse momentum (a), 
square of transverse momentum out of the event plane 
(b), alignment (c), and planarity (d). Estimation values 
for (p m ) and (p^) in case of 7r~p at 40 GeV/c were cal- 
culated on the basis of experimental data for eigenvalues 
of tensor S a ^ and for N c ^ from [23]. For 7r ± p interactions 
at initial momenta 4-25 GeV/c values of alignment and 
planarity corresponded to definitions used in present work 
(Tabled]) have been estimated on the basis of experimental 
results from [54]. 

One can see that behavior of (p? n )(N) (FigH^) and 
(Pg t )(AT) (Fig. [Sb) depends both on initial energy and 
from reaction type. Available experimental results for com- 
ponents of transverse momentum under study do not allow 
to make conclusions concerning presence and character of 
the universal analytical dependence describing experimen- 
tal data for different interactions and / or initial energies. 
But one can conclude that hadronic interactions show the 
decreasing of (pf n )(N) and increasing of (p^ ut )(N) with 
multiplicity increasing in general. 

Parameters characterized of alignment (Fig(SJ:) and 
planarity (FigJBJi) of event show smooth growth at mul- 
tiplicity increasing for all types of interactions and initial 
energy domain under study, i.e. the events are less planar 
and are characterized smaller alignment at larger hadron 
multiplicity in final state. Decreasing of values of param- 
eter (A) is observed with increasing of beam momentum, 
i.e. the alignment of hadrons along the principal axis ri\ 
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• Tr Ne, 6.2 GeV/c O K'p, 8.25 GeV/c a PP, 147 GeV/c ■ (v + v)p, 4.0 GeV 

T pp, 22.4 GeV/c <j. K'p, 32 GeV/c <> e + e", 57.8 GeV * (v + v)p, 5.0 GeV 

A 7t "P, 40 GeV/c ■ K + p, 70 GeV/c • (v + v)p, 2.8 GeV + (v + v)p, 6.5 GeV 

& 7t + p, 147GeV/c □ K + p, 147GeV/cx (V+v)p, 3.3 GeV * (V+v)p, 10.5 GeV 



Fig. 5. Dependence of average values of sphericity (a) and thrust (b) on multiplicity of identified particles in event. Experimental 
results for yr'Ne are from [70, 71], for K'p - from [22, 69], for 7r"p at 40 GeV/c - from [23], for K+p at 70 GeV/c - from [56], for 
pp at 22.4 GeV/c - from [67], for h + p at 147 GeV/c - from [63], for (u + v)p - from [48]. Estimations for an average thrust values 
in case of e+e~ are obtained on the basis of results from [90] . For (p + v)p average invariant mass of hadronic final state, (W) , is 
indicated above, for e + e~ - initial energy. Dotted thin lines marked as "1" and "2" at (a) are phenomenological approximations 
from [15]. Fit results from other publications, namely, for sphericity (a) and thrust (b) in h + p hadronic interactions from [63]; for 
e + e~ annihilation for thrust (b) from [90], are shown by thick dotted lines. Statistical errors are shown. See text for explanation 
of fit curves in more detail. 
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• 7i "Ne, 6.2 GeV/c X K + p, 32 GeV/c tc "P, 16.2 GeV/c £ 7t + p, 4 GeV/c 

* K " p, 8.25 GeV/c O p~p, 22.4 GeV/c * ji _ p, 25 GeV/c ■ tc + p, 8 GeV/c 
T K " p, 32 GeV/c A h 1 Be, 200 GeV/c A 7i"p, 40 GeV/c □ n + p, 1 6 GeV/c 



Fig. 6. Average values of collective parameters describing the event shape in the "principal axes" coordinate system depend on 
multiplicity of the secondary identified particles: a - (pf n ), b - (p ou t), c ~~ (A), d - (P). Values of (Pj n ) and (p ou ^) are shown 
in (GeV/c) 2 . Results for vr^Ne are from [70, 71], for 7r"p at 16, 25 GeV/c and for vr+p at 4, 8, 16 GeV/c are from [54], for K~p 
- from [22, 69], for pp at 22.4 GeV/c - from [67], for K+p at 32 GeV/c - from [68], for vr^p at 40 GeV/c - from [23], for h ± Be 
at 200 GeV/c - from [73]. Statistical errors are shown. See text for explanation of the curves. 
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Table 9. Fit results for multiplicity dependence in hadronic interactions 



Sample 




Fit parameter 








ai 


«2 0,3 


(fixed) 


X 2 /ndf 




7T _ NC 


0.97 ±0.05 


-1.44 ±0.13 


-0.5 


1.89 


PP 


0.91 ±0.07 


-1.9± 0.2 


-II- 


0.41 


7T~p 


0.680 ± 0.009 


-1.07 ±0.03 


-II- 


6.33 


K+p 


0.021 ±0.009 


0.028 ± 0.001 


1.0 


0.63 


(T) 


7T p 


0.594 ±0.015 


0.57 ± 0.04 


-0.5 


0.22 


K+p 


0.949 ± 0.005 


-0.015 ±0.001 


1.0 


1.94 



(jet axis) is amplified for higher initial energies. Influence 
of strange particles in a final state leads to small values 
of alignment in case of K~p interactions at initial mo- 
mentum 8.25 GeV/c [22] at small multiplicities. Moreover 
estimations of alignment parameter for K+p at 32 GeV/c 
[68, 69] are smaller significantly in the case of spheric- 
ity tensor definition Q than that for other experiments 
at close energies. Comparison of FigJBJ; and Fig(SJi demon- 
strates, that the parameter (P) is less sensitive to presence 
of strange particles in final state, than alignment. Detail 
additional study shows that behaviors of dependencies at 
Fig[5]and at Fig|6]are similar for total and charged multi- 
plicities and the general picture changes insignificantly at 
substitution TV^ot on -^ch an< ^ vase versa - 

The situation with theoretical description of behavior 
of multiplicity dependence for collective parameters differs 
from the situation for energy dependence of some observ- 
ables, at least. Apparently, there is no theoretical or phe- 
nomenological predictions for multiplicity dependence be- 
yond [15] or experimental data fits for main collective vari- 
ables [63, 90]. Taking into account monotonic and smooth 
behavior of dependencies under study, experimental data 
for average values of collective parameters in hadronic 
and (anti)ncutrino interactions have been fitted by gen- 
eral power function (JT]) for Y = S, T, p? Pout' ^, P 
and X = N in the case considered in this section. De- 
tail study of fit results for various data ensembles for all 
collective observables has shown, that better fit quality is 
achieved, usually, for a special case of (JT]) at the follow- 
ing fixed values = —0.5 or a 3 = 1.0. The first case 
agrees with function suggested in [15] for sphericity. Ap- 
proximations are shown by thick solid line for 7r _ Ne and 
by thin dashed line for pp at 22.4 GeV/c on FigJS^; by 
thick dashed line for 7r~p and by thin solid line for K+p 
at 70 GeV/c on FigEl Results for fit by function © at 
fixed a3 are shown in Table [5] for hadronic interactions 
under study. As seen, fit quality are acceptable for most 
data samples. One need to emphasize that general func- 
tion © allows to get the best quality for K+p collisions 
(x 2 /ndf = 0.01) with the following values of fit parame- 
ters: ai = 0.08 ± 0.03, a 2 = 0.009 ± 0.006, a 3 = 1.4 ± 0.3. 
But the 03 value is close to unit within errors and this 
general approximation agrees quite reasonable with linear 
one. As indicated above, modified both sphericity tensor 



and definition of corresponding observable from original 
papers which slightly differs from definition indicated in 
Table[T]are used for K~p interactions. Reducing of influ- 
ence of leading particles results in some faster intensifica- 
tion of spherical event shape with multiplicity increasing 
in comparison with behavior of (S)(N) dependence for 
other hadronic interactions at close energies for sphericity 
tensor definition Otherwise, results for K~p interac- 
tions at 8.25 GeV/c and 32 GeV/c agree with other ex- 
perimental points qualitatively but for the first case cor- 
responding fit quality is statistical unacceptable due to, 
possibly, small experimental errors. The following results 
have been obtained for fit of (S)(N) for K~p collisions at 
32 GeV/c by ©: a x = -0.94 ± 0.01, a 2 = 0.608 ± 0.009, 
a 3 = 0.346 ± 0.006 at x 2 / ndf = 1.33 (we do not show 
the curve at FigJS^, for more clearly picture). Multiplic- 
ity dependence of sphericity for (anti)neutrino beams are 
fitted by function with fixed 03 = —0.5 at all ener- 
gies with statistically reasonable qualities. Results of fit 
for (v ± v)p interactions are presented in Table ITOl (we do 
not show corresponding curves at Fig(5£i for more clearly 
picture). 

Experimental data for various interaction shown at 
Figj5] and fit results for main collective variables (Tables 
l9"land[TU| allow to suggest that intensification of spherical 
event shape with multiplicity increasing seems more pow- 
erful in energy domain y/s > 11 GeV than that at lower 
initial energies. 

Table 10. Fit results for multiplicity dependence of (S) at 
(anti)neutrino beams 



(W), GeV 


Fit parameter 




ai 


ai 


X 2 /ndf 


4.0 


1.09 ±0.11 


-1.7 ±0.3 


1.34 


5.0 


1.08 ±0.12 


-1.8 ±0.3 


0.42 


6.5 


0.69 ±0.07 


-1.0 ±0.2 


0.19 


10.5 


0.72 ±0.11 


-1.4 ±0.3 


1.22 



Notations for approximation curves at Fig|5] are the 
same as well as on Figj5]for corresponding reactions. More- 
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Table 11. Values of fit parameters for function Q 



Data 


Fit parameter 




sample 


ai 




x 2 /ndf 


(Pi) 


vr~Ne 
K-p 

7T p 


0.082 ± 0.002 

0.06 ±0.02 
0.065 ± 0.004 


0.0 (fixed) 
0.15 ±0.06 
0.123 ±0.012 


0.88 
1.47 
2.78 






> 




7r~Ne 
K-p 

7T p 


0.045 ± 0.004 
0.104 ±0.004 - 
0.067 ± 0.002 - 


-0.07 ±0.01 
-0.190 ±0.008 
-0.102 ±0.005 


1.36 
1.98 
2.82 


(A) 


7r Ne 
PP 

7T~p 


0.77 ±0.06 
0.65 ±0.01 
0.50 ±0.01 


-1.17±0.14 
-1.10 ±0.03 
-0.79 ±0.03 


0.89 
3.96 
6.15 


(P) 


7r"Ne 

A"-beams 

7r-beams 


0.68 ±0.06 
0.885 ± 0.013 
0.912 ± 0.010 


-0.94 ±0.15 
-1.56 ±0.03 
-1.64 ±0.02 


2.06 
2.47 
8.13 



over the fit curve for K~p at 32 GeV/c is shown by solid 
thin line on FigfBk, b; for K~p at 8.25 GeV/c - by dot- 
ted line at FigJBJ;; and curve for all available hadronic in- 
teractions is demonstrated by thin solid line at FigfBJi. 
For multiplicity dependence of alignment in K~p reac- 
tion at 8.25 GeV/c the following values of fit parameters 
have been obtained for ([7]): a\ = —3.035 ± 0.008, a 2 = 
2.676 ±0.008, a 3 = 0.113 ±0.001, and xV n df = 1.55. Fit 
result for various experimental data samples by specific 
case of function (Q at fixed 03 = —0.5 are shown in Table 
1111 for all interactions and event shape variables with the 
exception of (p? )(N) for 7r~Ne at 6.2 GeV/c discussed 
below. As seen, multiplicity dependence for all collective 
variables are fitted by specific case of ([7J with reasonable 
quality in the case of 7r~Ne hadron-nuclear collisions. For 
this reaction dependence of (p 2 n ) on N (Figj6^,) is fitted 
with best quality by the simplest case of ([7]) which corre- 
sponds to constant. On the other hand fit by functional 
form (UJ) with fixed 03 = —0.5 shows reasonable and only 
slightly poorer quality (x 2 /ndf=1.02). Thus the behavior 
of (pf n )(N) does not contradict with specific case of func- 
tion J7J indicated above for hadron-nuclear reaction under 
study. The reasonable quality has been obtained for fit of 
data sample for (p? ) in K~p at 32 GeV/c without first 
point only. Dependencies of (pj n ) and (p^ut) on ^ shown 
for K ~p at 8.25 GeV/c on Fig [5^ and Fig[5p, respectively, 
agree with approximation function flTJ at a 3 = —0.5 qual- 
itatively but the corresponding fit qualities are statisti- 
cal unacceptable for this reaction due to dispersion of 
experimental points and, possibly, small errors. As seen 
from Table 1111 planarity results obtained in interactions 
with kaon beams are fitted by (J7J with reasonable % 2 /ndf. 
Experimental results for pion beams agree with approx- 



imation function ([7]) qualitatively only and fit quality is 
some poorer for this incoming particle than that for kaon 
beams. The approximation curves for two beam types un- 
der study are very close to each other, moreover results 
for all hadronic interactions demonstrate similar values 
for planarity parameter (P) at corresponding multiplici- 
ties (Fig|6ji). Thus we have attempted to fit all available 
experimental points by {7J at fixed (13 — —0.5. One can 
conclude from FigJSJi that total sample of experimental 
results for hadronic interactions agrees with specific case 
of function (UJ at qualitative level but the fit quality is 
statistically unacceptable (x 2 /ndf ~ 20). 

Additional investigations show the following. First of 
all, power function ([7J and it's special case at fixed a 3 = 
—0.5 describe dependencies of collective variables on char- 
ged multiplicity reasonable at qualitative level too. Sec- 
ond, as expected, using of modified sphericity tensor with 
reduced momenta of final state particles for K~p [22, 69] 
result in growth of values of (A) (especially) and (P) with 
the preservation of functional form for corresponding de- 
pendencies. On the other hand the behavior of depen- 
dence of components of particle transverse momentum in 
and out of the event plane on multiplicity changes dra- 
matically, namely, (p 2 n ) and (p^ut) are snown the power 
growth with N increasing and significantly larger values 
than that for sphericity tensor definition 

6 Summary 

In the framework of present paper the database of ex- 
perimental results for traditional event shape variables 
was created for wide set of interaction types. This database 
is used for investigation of dependencies of collective pa- 
rameters on initial energy and multiplicity in various in- 
teractions and for joint samples. 

At intermediate energies influence of nuclear matter, 
in particular, multinucleon interactions leads to signifi- 
cantly more isotropic distribution of secondary particles 
for hadron-nuclear reactions in comparison with hadron- 
hadron ones, to growth of transverse momentum, attenu- 
ation of alignment and planarity of events. Leading par- 
ticles result in the inverse influence on event shape in 
hadronic reactions in comparison with intranuclear rescat- 
tering. 

Energy dependencies of the main collective observ- 
ables sphericity and thrust show universal behavior for 
various interactions. Phenomenological models based on 
the QCD, in particular, dispersive model allow to describe 
(T) (s) in wide range of initial energies down to strongly 
non-perturbative domain y/s ^2 — 3 GeV at qualitative 
level at least. Approximation of energy dependence of av- 
erage thrust values for various interactions allows to ex- 
tract as(Mz) estimations which are in agreement both 
with world average value and with results obtained by 
other methods. Values of transverse momentum compo- 
nents which are in / out of event plane for e + e~ interac- 
tions increase faster with growth of ^fs, than that for ep 
and hadronic interactions. Energy dependencies of align- 
ment and planarity suppose the descriptions by universal 
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empirical functions for different interaction types. Using 
suggested analytic approximation functions estimations of 
values of collective parameters under study have been ob- 
tained for present and future collider-research facilities. At 
y/s ^5 — 7 TeV average values of collective variables stud- 
ied in the present paper do not depend on initial energy 
within errors. Thus, these values of y/s can be considered 
as an estimation of the low bound of asymptotic region 
for traditional collective parameters. 

Spherical event shape is amplified with multiplicity in- 
creasing and one can suggest, that the universal estima- 
tion of the low energy boundary for experimental appear- 
ance of event jet structure in multiparticle production pro- 
cesses is y/s^ ~ 3 GeV. Usually, multiplicity dependence 
of collective variables under study agree with power func- 
tion in energy domain y/s < 12 GeV at qualitative level 
at least. 
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